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Abstract The origin of the extragalactic gamma-ray background (EGRB) is still an open 
question, even after nearly forty years of its discovery. The emission could originate from 
either truly diffuse processes or from unresolved point sources. Although the majority of 
the 271 point sources detected by EGRET (Energetic Gamma Ray Experiment Telescope) 
are unidentified, of the identified sources, blazars are the dominant candidates. Therefore, 
unresolved blazars may be considered the main contributor to the EGRB, and many stud- 
ies have been carried out to understand their distribution, evolution and contribution to 
the EGRB. Considering that 7-ray emission comes mostly from jets of blazars and that 
the jet emission decreases rapidly with increasing jet to line-of-sight angle, it is not sur- 
prising that EGRET was not able to detect many large inclination angle active galactic 
nuclei (AGNs). Though Fermi could only detect a few large inclination angle AGNs in 
the first three months' survey, it is expected to detect many such sources in the near future. 
Since non-blazar AGNs are expected to have higher density as compared to blazars, these 
could also contribute significantly to the EGRB. In this paper we discuss contributions 
from unresolved discrete sources including normal galaxies, starburst galaxies, blazars 
and off-axis AGNs to the EGRB. 
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1 INTRODUCTION 

Extragalactic diffuse emission has been observed nearly at all wavelengths ranging from radio to 
gamma-rays. Nearly 18 years ago, gamma-ray astronomy received a substantial boost with the launch 
of the Compton Gamma-Ray Observatory (CGRO) in 1991 (Kanbach et al. 119881 ). With significantly 
improved sensitivity over previous experiments and a low instrumental background, EGRET (on board 
CGRO) provided a new platform to study the spectrum and distribution of the extragalactic emission in 
greater detail than it was possible in the past. In 2008, with the launch of the Fermi Gamma-ray Space 
Telescope, one expects another significant leap in determining the nature of the extragalactic gamma-ray 
background (EGRB), and first results from Fermi have already been quite remarkable (Abdo et al. 120091 1. 

To date, the knowledge of the EGRB intensity and spectrum (Sreekumar et al. 119981 Strong, 
Moskalenko & Reimer 120041 ) is derived from the all-sky data from EGRET spanning more than four 
years of observations. This will soon be updated based on new results from Fermi all-sky observations. 
Sreekumar et al. i 19981 ) used the EGRET data to derive the extragalactic emission as a constant com- 
ponent of the total observed emission that is uncorrelated with line-of-sight column density of matter 
for thirty-six independent regions of the sky. The average spectrum of the EGRB was found to be well 
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fit with a single power-law characterized by a spectral index of (-2.10 ± 0.03) and the integrated flux 
above 100 MeV was derived to be (1.45 ± 0.05) x 10~ 5 photons cm~ 2 s _1 sr _1 . In a different anal- 
ysis, Strong, Moskalenko & Reimer (2004) used a more extensive cosmic -ray model incorporated in 
the GALPROP code, inferred a new model for Galactic diffuse continuum gamma-rays, and found the 
EGRB integrated flux above 100 MeV to be (1.11 ± 0.01) x 10 -5 photons cm~ 2 s _1 sr~\ slightly 
lower than that obtained by Sreekumar et al. (119981 1. In this case, the EGRB spectrum if fitted with a 
single power law, was found to yield an index of -2.17 ±0.04 with a possible spectral break around 1.5 
GeV. 

The origin of the extragalactic emission has long been a subject of study (see Kneiske (|2008|> for a 
recent review). As in the case of diffuse emission at most wavelength bands, it can be interpreted either 
as truly diffuse emission of cosmological origin or as arising from limited capability of instruments to 
detect large number of weak point sources. Truly diffuse emission can arise from numerous processes 
such as black hole evaporation, particle acceleration by intergalactic shocks produced during large scale 
structure formation (Loeb & Waxman |2000| Gabici & Blasi 2003), for example. In the alternate scenario 
wherein unresolved sources contribute, the mean intensity of diffuse emission should decrease with in- 
creasing angular resolution of instruments. The final EGRET point source catalog (Hartman et al. 11999b 
contains 271 sources (>100 MeV), which includes five pulsars, one probable radio galaxy (Cen A), 66 
high confidence identifications of a sub-class of active galactic nuclei called blazars and a single external 
normal galaxy, the Large Magellanic Cloud (LMC). In addition, 27 low confidence potential blazar iden- 
tifications are listed. Therefore, blazars could contribute significantly to the EGRB. Studies published 
so far estimate that the blazar contribution to EGRB can range from 25% to 100% (Stecker, Salamon 
& M alkan [19931 Chiang et al. [19951 C hiang & Mukherjee [19981 Stecker & Salamon [19961 Mwcke & 
Pohl[2000, Narumoto & Totani. 12006b . Similar contributions from other source classes have also been 
reported (Stecker 120071 Thompson et al. 120071 Pavlidou et al. 120071 Bhattacharya & Sreekumar |2009l . 

EGRET observations have provided strong arguments for AGNs as a source class that can contribute 
to the unresolved point source component of EGRB; however, uncertainties in the current predictions 
provide justifcation to consider contributions from other gamma-ray bright sources. In this paper, we 
re-examine the contribution from blazars to the EGRB using updated source catalogs (see details in 
Bhattacharya, Sreekumar & Mukherjee (120091 1). We also include an estimate of the contribution from 
off-axis AGNs to the EGRB. While we await new results from Fermi on the EGRB, the analysis pre- 
sented here for the EGRET data will be an useful exercise for future such studies to be carried out with 
the Fermi data. 



2 RADIO-LOUD UNIFICATION SCHEME & OFF-AXIS AGN 

AGN studies in the past have led to various source classifications. There have been systematic efforts 
to explain the observed differences amongst the various classes of AGNs, including the more popular 
schematic models known as "Unification schemes" (see e.g. Antonucci (119931 1; Urry & Padovani d 19951 ); 
Robson ( 1996 )). A central feature of current unified schemes is that observed properties, and thus classi- 
fication of a given AGN depend on its orientation. According to the Fanaroff & Riley (FR) classification 
(Fanaroff & Riley |1974l > of AGN, radio loud quasars have FR II morphology (radio emission brighter in 
the lobes or outer parts of the jets than in the central core), while BL Lac objects have FR I morphology 
(radio emission peaks in the central core of the galaxy). It is assumed that steep-spectrum and flat- 
spectrum radio quasars (SSRQs and FSRQs) and FR II radio galaxies are from one parent population, 
while BL Lacs and FR I radio galaxies are from another parent population. Further, it is assumed that for 
FSRQs, the line-of-sight is nearly aligned to the jet axis. Therefore, strong relativistic Doppler beam- 
ing of the jet emission produces highly variable and continuum-dominated emission. For SSRQs, the 
line-of-sight is slightly outside the jet cone; however, the continuum remains beamed albeit weakened. 
As one moves away from the jet axis, the central continuum flux falls and the nucleus appears more 
like a FR II BLRG (broad-line radio galaxy). If the dusty torus obscures the line-of-sight, the source 
looks like a FR II narrow-line radio galaxy (NLRG). Similarly, for BL Lacs, the jet is nearly aligned to 
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the line-of-sight and strong relativistic Doppler beaming is expected. Away from the jet axis, one will 
observe FR I type galaxies. 

According to the above unification scenario, AGNs having very small jet to line-of-sight angle are 
termed "blazars." When the AGN is viewed at small angles to the jet axis, one expects greater probability 
to detect high energy (X-ray and gamma-ray) emission. Almost all the AGNs detected by EGRET 
above 100 MeV are classified as blazars. It is assumed that the high energy flux from a jet decreases 
very rapidly with increasing jet to line-of-sight angle. So most AGNs with large jet inclination angles 
are expected to be characterised by small gamma-ray flux and hence may not have been detected by 
EGRET. However, EGRET detected the radio galaxy, Cen A (~ 70° inclination angle; Graham ( 11979b ; 
Tingay et al. d!9981 >; Hartman et al. ( 119991 and possibly, NGC 6251 (~ 45° inclination angle; Sudou 
& Taniguchi (2000); Mukherjee et al. ( 120021 )), which have large inclination angles. Although the Fermi 
Bright Source Catalog lists only one more radio galaxy, NGC 1275, in 7-rays (Abdo et al. |2009j, it is 
expected to detect many such radio galaxies in the near future. Considering the unification model, the 
number of AGNs with increasing inclination angles is expected to be much higher than the number of 
blazars. Though the high-energy jet emission is expected to decrease rapidly with increasing angles, 
considering the larger number densities over blazars, these sources could contribute significantly to the 
EGRB. The contributions of dim blazars and radio/gamma-ray galaxies to the EGRB has recently been 
studied by Dermer d2007bl l. Here, we present a calculation of the contributions from radio-loud quasars 
and FR I and FR II galaxies (relative to blazars) to the EGRB. 

3 EMISSION FROM AGN JET 

Gamma-ray emission in blazars is believed to originate in relativistically beamed jets, thought to be 
powered by accretion of matter onto a central supermassive black hole. The AGN jet is assumed to 
be made up of a large number of discrete blobs, consisting of energetic particles that may be either 
electrons or hadrons ("leptonic" vs "hadronic" blazar models; see for e.g. Boettcher (2007); Reimer, 
Joshi & Bottcher (2008 ) for a review). It is assumed that as the blob propagates along the jet it slows 
down. In leptonic models, for example, the emission is mainly due to the synchrotron losses of the 
electrons and the inverse Compton scattering of either these synchrotron photons from within the jet 
("Synchrotron Self-Compton" or SSC model; Maraschi et al. (11992b ; Bloom & Marscher ( 119961 1, for 
example) or the accretion disk photons ("External Compton" or EC model; see e.g, Dermer, Schlickeiser 
& Mastichiadis ( H992K hereafter DSM). Several other alternative models exist, and the exact nature of 
the gamma-ray emission mechanism in blazars is still unresolved. 
We model the gamma-ray luminosity L of an AGN at an inclination angle 9 as 



where, L max is the maximum luminosity, and gives the angular dependence of the luminosity, with 
£,max = 1- The on-axis luminosity (Lq) is given by L max x £(0). 

We consider the following two broad categories of models in the subsequent sections: 

a) Synchrotron Self Compton emission (SSC) model. 

b) External Compton scattering of accretion disk photon (EC) model, 
a) Synchrotron Self Compton emission (SSC) model: 

Let T be the bulk Lorentz factor of the blob. As it traverses the jet, V varies from its initial value (T max ) 
to some minimum value (r m ;„) below which high-energy emission ceases. It is considered that the 
decay length scale is much smaller than the source luminosity distance. For the source with a power law 
spectrum F{y) oc v a the transformation can be written as (Kembhavi & Narlikar [l999l ) 



L{6) = L. 



'max 



x£(0) 



(i) 



F{y) = ( 
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l + z 
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Here F{v) is the flux at a frequency v, a is the spectral index, and z the redshift of the source. 
The primed quantities are in the blob rest frame and unprimed quantities are in the observer frame. The 
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Doppler boost factor S is expressed as 

s = r(i - pcose) (3) 

where is the angle between the jet axis and the line-of-sight, and (3 is the speed of the blob (in the unit 
of velocity of light). 

It is considered that the power law energy distribution of electrons in the blob is given by n e (7) = no7~ p 
for 1 < < 71 < 7 < 72 , and zero otherwise. Here, 7 is the Lorentz factor of the electrons. If the electrons 
are emitting synchrotron radiation with a spectral index a and electron energy index p, then a and p are 
related by a = For the sample of EGRET-detected gamma-ray blazars, Bhattacharya, Sreekumar 
& Mukherjee d20Q9b found that the average a value is 1.34 for FSRQs and 1.08 for BL Lacs. The 
corresponding p values are 3.68 and 3.16 for FSRQs and BL Lacs, respectively. The angle-dependent 
part of luminosity is proportional to 6^ (Dermer [l9951 l. 
b) External Compton (EC) Model 

In a different scenario, accretion disk thermal photons can enter into the blob and can be upscattered by 
high energy electrons (DSM 1992). The following assumptions are made: 

1 . The relativistic electrons are homogeneously distributed throughout the blob. 

2. The randomly-oriented magnetic field has uniform strength throughout the blob. 

3. In the blob frame (BF), electron energy distribution is isotropic. 

4. The accretion disk and the core of the AGN, emit photons which are isotropic in the observer frame 
(OF) ( photon flux oc \, where, r is the distance of the blob from the central source ) . These photons 
enter the outflowing blob from behind (upstream). 

5. Photons pass through the blob following trajectories parallel to the jet axis. This assumption is justified 
if the bob is sufficiently far away from the central source. 

6. The blob is optically thin to Thompson scattering along the jet axis. This assumption is necessary in 
order to produce highly polarized radio emission (DSM (1992) and references therein). 

The angle-dependent part in the flux and hence in the luminosity is proportional to <5 3+p (l — /x^) 2 ^ - , 
where // s = cosO, 9 is the jet to line-of-sight angle (DSM 1992; Weferling & Schlickeiser ll9991 l. 

3.1 Angular dependence (£(6>)) of the luminosity 

In order to incorporate the slowing down of the jet, the empirical relation given by Georganopoulos and 
Kazanas j2003aH2003bl ) is used. 

r(r) = T max if < r < r 

ill- I'll (4) 

T(r) = T x ( — y 2 if r min <r< r max (5) 

Here, r is the distance of the blob from the central source, and from X-ray data it can be shown that 
Tq = 3 for ro = 1 kpc (Georganopoulos and Kazanas !2003a). Since any source no matter how faint it 
is, will contribute to the background, T m i n is assumed to be 1 . The values of r m i n and r max correspond 
to Y max and T mm , respectively. 

Hence, for the SSC model, the integrated luminosity in the observer frame (OF) is given by 

^ L s s F c dr = r maX L% S F C xS^dr 
Jo 

= K ssc x I ssc {0) (6) 

Similarly, for the EC model, 



L E c F dr = 



EC 
J BF 



(J) 
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The I ssc (9) and I EC (0) give the angle dependence of the luminosity while K ssc and K EC are 
^-independent. Since we are mainly interested in how the jet emission changes with angle, we only 
focus on I ssc (8) and I EC (6). The function £(9) has been defined before in such a way that for the 
SSC model, 



max 

and for the EC model 



e C {0) = (9) 



e°{0) = do) 

-^inax 




20 40 60 80 20 40 60 80 

Jet to line-of-sight angle (8) Jet to line-of-sight angle (8) 



Fig. 1 for different emission scenarios 

(discrete jet model) forFSRQs. The solid lines 
correspond to situations where plasma blob 
decelerates and hence T changes. T max = 10 
and T min = 1. Thick line = EC model (DSM 
1992) and thin line = SSC model (Dermer 
1995). The dashed lines correspond to V = 
constant = 10. The vertical line corresponds to 
10° viewing angle. See the electronic edition 
of the Journal for a color version of this figure. 



Fig. 2 £(6) for different emission scenarios 
(continuous jet model) for FSRQs. The solid 
lines correspond to situations where plasma 
blob decelerates and hence V changes. T max 
= 10 and T m i n = 1. Thick line = EC model 
and thin line = SSC model. The dashed lines 
correspond to T = constant =10. The verti- 
cal line corresponds to 10° viewing angle. See 
the electronic edition of the Journal for a color 
version of this figure. 



4 CONTRIBUTION TO THE EGRB FROM AGNS 
4.1 Blazar contribution 

We begin with the contribution to the EGRB from blazars. Previous works on this can be broadly clas- 
sified into two groups. 
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Fig. 3 £(0) for different emission scenarios 
(discrete jet model) for BL Lacs. The solid 
lines correspond to situations where plasma 
blob decelerates and hence T changes. Y max = 
4 and Y mm = 1 . Thick line = EC model (DSM 
1992) and thin line = SSC model (Dermer 
1995). The dashed lines correspond to T = 
constant = 4. The vertical line corresponds to 
10° viewing angle. See the electronic edition 
of the Journal for a color version of this figure. 



Fig. 4 £(0) for different emission scenarios 
(continuous jet model) for BL Lacs. The solid 
lines correspond to situations where plasma 
blob decelerates and hence Y changes. Y max 
= 4 and Y m i n = 1. Thick line = EC model 
and thin line = SSC model. The dashed lines 
correspond to Y = constant = 4. The vertical 
line corresponds to 10° viewing angle. See the 
electronic edition of the Journal for a color 
version of this figure. 



(a) A linear relationship is assumed between the gamma-ray luminosity of blazars and the luminosity 
at some other wavelength so that the gamma-ray luminosity function can be scaled with the luminosity 
function of that wavelength (Stecker, Salamon & Malkan 1993; Stecker & Salamon ll996l Narumoto & 
Totani.|2006l>. 

(b) Constructing the gamma-ray luminosity function of blazars using the gamma-ray catalog (Chiang et 
al. [19951 Chiang & Mukherjee [1998|>- 

The flux Fo (per solid angle) coming from on-axis sources can be expressed as 

4tt7 dz J Lmtn0 AttD l 2 

where Lq and <f>o(Lo, z) are the on-axis luminosity and on-axis luminosity function, respectively. 

Stecker, Salamon & Malkan (1993) and Stecker & Salamon (|19961 l adopted the first approach and 
assumed a linear relationship between the radio and the gamma-ray luminosity of blazars. They found 
that unresolved blazars can fully explain the observed level of EGRB. In a similar approach, Narumoto 
& Totani (2006) considered a linear relationship between the X-ray and the gamma-ray luminosity of 
blazars and found unresolved blazars contribution to the EGRB to be 25-50%. Using approach (b), 
Chiang et al. d 19951 > and Chiang & Mukherjee ( 11998t constructed the gamma-ray luminosity function 
from EGRET-detected blazars alone and found the contribution to be ~ 20 %, using a catalog of 34 
sources. 

Recently, Bhattacharya, Sreekumar & Mukherjee (|20091 > re-calculated the luminosity functions of 
FSRQs and BL Lacs and their evolution, but using the most recent EGRET gamma-ray source identifi- 
cation catalog of (Sowards-Emmerd, Romani & Michelson 2003; Sowards-Emmerd et al. 120041 1 which 
almost doubled the blazar count from the EGRET 3EG catalog, and considering FSRQs and BL Lacs as 
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separate classes of objects. In that study, BL Lacs were found not to show any evidence for evolution. 
The luminosity function was found to be best described by a single power law with index -2.37 ± 0.03, 
and it was not possible to determine if there was any break in the luminosity function. For FSRQs, the 
< v v > test indicated strong evolution, and a pure luminosity evolution was considered. Both ex- 
ponential and power law evolution functions were examined. A broken power law form of de-evolved 
luminosity function was assumed (lower luminosity end index oti, upper luminosity end index 0:2 and a 
break luminosity L b)- Table 2 of Bhattacharya, Sreekumar & Mukherjee ( 120091 1 shows values of lumi- 
nosity function parameters for two different evolutionary models. The source distribution with redshift 
is better explained for an exponential evolution function than a power law function. For the exponential 
evolution function, the contribution to EGRB becomes ~9 xl0~ 7 ph cm 2 s _1 which is ~ 6% of the 
total EGRB as calculated from the EGRET data. The maximum contribution (considering error limits 
from Table 2 of Bhattacharya, Sreekumar & Mukherjee ( 120091 )) is estimated as ^15%. 

For BL Lacs we incorporated the break luminosity and lower luminosity index from that of FSRQs. 
The calculation constrains the EGRB contribution to < 6 %. 

Early results from Fermi also found no evolution for BL Lacs (Abdo et al. 120091 ), though their 
BL Lac luminosity function is found to be harder than that calculated by Bhattacharya, Sreekumar 
& Mukherjee (2 009 ) from the EGRET data. There may be several factors behind these differences. 
Since these sources are strongly time-variable, the time-averaged source characteristics are needed to be 
considered in order to address contributions of source classes to the background. This should be possible 
from the next few years of Fermi data. 



4.2 Off-axis AGN contribution 

The on-axis source contribution to the EGRB can be estimated from Eq. QT| For a particular angle 8, 
L(8) is linearly related to Lq (Eq.Q}, so the luminosity function of off-axis sources can be replaced by 
the on-axis luminosity function. Hence the contribution to the EGRB from AGNs having jet inclination 
angle 8, can be expressed as 



F(8) = ±- r a *^dz [ LHmMZ MLe,z) Le{ ] + *V *- dL ° 
4tt J dz J Lmin>e At:D l z 

4tt J dz J Lmin , f (0) x 47rD L 











——dz 






dz 












dz 






dz 


'L min , 


1 






-in 


r 

Jo 


dz dZ j L 



<p (L ,z) 2 dLo (14) 

So the total contribution to the EGRB from AGNs inclined at any angle, 

F = C° F{8)d8 
that is F = jkpr Jq® £(8)d8x contribution from on-axis sources. 

Here, it is considered that EGRET detected blazars have jet inclination angle between to 10°. So the 
contribution from the on-axis sources as detected by EGRET is 

F EG ret= / F{6)d6 (15) 
and the contribution from the off-axis sources will be given by 



^.90° 

Fqff = / F{8)d8 
Jio° 



(16) 
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Therefore, 



Fqff = C W)dO 
Fegret Cmde 



and the off-axis AGNs contribution to the EGRB is 



90° 

Fqff = J A°o° x Fegret (18) 



Similarly, considering a jet to line-of-sight angle (8) range from 90° to 180°, one can include counter 
jet contribution also, although this is much smaller than the jet contribution. 

We examined the angular dependence of jet emission in AGNs for two different emission models for the 
discrete jet model as discussed above. The situation of a continuous jet(see e.g., Urry & Padovani (119951 1) 

is also examined with the angle dependence only coming through the Doppler beaming factor. For EC 

p+i 

emission of the DSM model, the peak emission is shifted from zero angle due to the (1 — /i s ) 2 factor. 
For a continuous jet and EC emission model, and considering the angle dependence of the luminosity to 
be only due to the Doppler beaming factor, the angle dependent part of the luminosity is proportional to 
5 2+p . However for SSC emission model, the angle dependent part of the luminosity is proportional to 

S 2+a . 

Our results are presented in figures 1-4. Fig.[T]andFig.[2]show the angle dependence of the luminosity for 
FSRQs for the discrete jet model and for the continuous jet model, respectively. For FSRQs we consider 
Fmax to be 10 (Jorstad et al. (120011 ); Dermer (2007a) and references therein). It is found that for external 
Compton model (from DSM model) the off-axis contribution is ~ 15% of the FSRQ contribution. We 
give our results in Table 1 . 

Fig. [3] and Fig. |4]show the angle dependence of the luminosity for BL Lacs for the discrete jet model 
and for the continuous jet model, respectively. For BL Lacs, T max is assumed to be 4 (Dermer J2007a| ) 
and references therein). Since the T max is less than that for FSRQs, the jet emission is less beamed, 
and hence the off-axis contribution (relative to on-axis) is larger. If a continuous jet is considered, then 
the off-axis contribution can be ~ 29% of the BL Lac contribution. For the DSM model and the EC 
emission model, the off-axis contribution is larger than the on-axis contribution (TableQ]). 
We discuss our results in the context of previous work in this area. Miicke & Pohl ( 120001 also estimated 
the contribution from FSRQs and BL Lacs separately. Considering a leptonic jet model, they calculated 
the non-thermal emission from blazars, with the inverse Compton scattering of soft photons from the 
accretion disk (EC process) dominating the synchrotron self Compton process. They estimated the con- 
tribution from these source classes to the EGRB in the context of the AGN unification paradigm, and 
found the total contribution from unresolved blazars (both "aligned" and "misaligned") to the EGRB as 
(20-40)%. In a similar approach, Dermer d2007al ) used a physical model to fit the redshift and size distri- 
bution of EGRET detected blazars. Dermer d2007al ) estimated that the FSRQs and BL Lacs contribution 
to the EGRB at 1 GeV as ~(10-15)% and ~(2-3)%, respectively. The approach adopted by Mitcke & 
Pohl (2000) and Dermer (2007a) considered both "aligned" and"misaligned" blazars. Cao & Bai ( 120081 ) 
calculated the contribution from RLQs and FR II galaxies to the EGRB in a different approach. Cao 
& Bai (120081) derived the parent radio luminosity function of radio quasars/FR lis from the FSRQ lu- 
minosity function derived by Padovani et al. J20071 I. They assumed that gamma-rays in AGN jets are 
produced by SSC and EC processes. They further assumed a power law distribution of the bulk Lorentz 
factors of jets. Using the derived parent luminosity function, they calculated the observed gamma-ray 
luminosity functions for FSRQs, SSRQs and FR II galaxies, and their contribution to the EGRB. Cao & 
Bai (f2008) estimated the contribution from all radio quasars/FR II galaxies to be ^30% of the EGRB. 
They did not calculate the contribution from BL Lacs/FR I galaxies to the EGRB. Stawarz, Kneiske 
& Kataoka d2006t calculated the gamma-ray emission from FR I galaxies assuming that the observed 
X-ray emission in the bright knots of the kiloparsec scale jets are due to synchrotron emission. Including 
the contribution from secondary gamma-ray photons, they found that FR I galaxies contribute ~1% to 
the EGRB. 
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Table 1 Off-axis contribution to the EGRB (relative to on-axis) of FSRQs & BL Lacs for 
different jet models. 



Source Type Emission Models discrete jet Continuous Jet 
FSRQ EC Model ~ 15 % [T, 1] ~ 0.1 % [T, 3] 

(A max — 10) 

a = 1.34 

SSC Model ~ 1 % [T, 2] ~ 3 % [T, 3] 



BL Lac EC Model > 100 % [T, 1] ~ 8 % [T, 3] 

(J-^max — 4) 

a = 1.08 

SSC Model ~ 15 % [T, 2] ~ 29 % [T, 3] 



References: (T) This paper; (1) Dermer, Schlickeiser & Mastichiadis (1992); (2) Dermer (1995); (3) Urry & 
Padovani (1995) 

As noted earlier (Section 4.1), we find that the blazar contribution can go upto ~ 21% of EGRB. If 
one considers that the contribution from non-blazar radio-loud AGNs to the EGRB is ^(30-40)% of 
the blazar contribution, the maximum contribution from all radio quasars, FR I and FR II galaxies 
becomes ~ 29% of EGRB. With better sensitivity and resolution, Fermi is expected to detect more 
AGNs (~ 1000; Dermer d2007al ): Cao & Bai ( 120081 )). which will provide a better estimation of the 
contribution from those AGNs to the EGRB. The substantial increase in detected sources within various 
source classes could also permit a more direct determination of the luminosity function and evolution of 
off-axis AGNs. 

5 RESULT & DISCUSSION 

In an earlier paper (Bhattacharya & Sreekumar |2009l ). we found that the contribution from normal galax- 
ies is ~ 1% and that from starburst galaxies is ~5% (maximum ~6%) of the total EGRB. In this paper, 
we found that the contribution from unresolved AGNs is ~15% of the EGRB. So, we conclude that the 
total contribution from normal galaxies, starburst galaxies and blazars to the extragalactic gamma-ray 
background is ~ 21%. If the upper limit to the non-blazar radio-loud AGNs contribution is ^40% of the 
blazar contribution, the total contribution from normal galaxies, starburst galaxies and radio-loud AGNs 
to the EGRB can atmost be ~36%. Truly diffuse emission processes can at present account for only a 
smaller fraction of EGRB. 

Most of the identified EGRET sources are FSRQs. Their average spectra have a spectral index of 
2.34±0.15, where the EGRB is well described by a single power law of index 2.10±0.03 (Sreekumar 
et al. 119981 ). This spectral incompatibility suggests the need for source classes with flatter spectra. BL 
Lacs have an average spectral index of 2.08. But their luminosity function is not well determined due to 
limited number of detection of these sources by EGRET. 

Pavlidou et al. ( 2007 ) found that the contribution from the unresolved gamma-ray sources of the same 
class of unidentified EGRET sources, contribute significantly to the EGRB. Further improvement is 
only possible through the deeper understanding of nature and distribution of these source classes which 
is possible with Fermi. 

The new sky survey from Fermi is expected to give a much improved picture of EGRB. From this work 
it is evident that the unresolved sources of currently known gamma-ray source classes are not sufficient 
to explain fully the observed EGRB. Following the same methodology with Fermi data one can have 
better constrained luminosity functions of these source classes, construct luminosity functions of newly 
discovered source classes, address if there remains a spectral incompatibility with sources and revisit 
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the most interesting question of all " Is there a truly diffuse gamma-ray background component over 
and above the unresolved source contributions ?" 
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